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Abstract. - We compute the electron spin resonance (ESR) intensity for one-dimensional 
quantum wires in semiconductor heterostructures, taking into account electron-electron interac- 
tions and spin-orbit coupling. The ESR spectrum is shown to be very sensitive to interactions. 
While in the absence of interactions, the spectrum is a flat band, characteristic threshold sin- 
gularities appear in the interacting limit. This suggests the practical use of ESR to reveal spin 
dynamics in a Luttinger liquid. 



The search for the Luttinger liquid (LL) behavior of one-dimensional (ID) interacting 
fermions Jl]-|i|] has a rather long history. Recent experimental work has reported strong evi- 
dence via the tunneling density of states in single-wall carbon nanotubes || [|, for the edge 
states in the fractional quantum Hall effect Q] (which supposedly are chiral LLs not further 
discussed here), and via the conductance H and resonant tunneling in a quantum wire 
(QW) in semiconductor heterostructures. These studies have so far only probed the charge 
dynamics and are intrinsically insensitive to spin dynamics. However, spin dynamics can re- 
veal some of the most interesting and spectacular aspects of the non-Fermi liquid nature of a 
LL, in particular the phenomena of spin-charge separation and electron fractionalization 
In this Letter, we suggest to probe the spin dynamics of a LL via electron spin resonance 
(ESR), and provide detailed theoretical predictions for the ESR spectrum in an interacting 
ID QW. ESR spectra have already been obtained experimentally for the 2D electron gas in 
heterostructures |jl(|, and new fabrication advances such as cleaved-edge overgrowth tech- 
niques [jll| indicate that such experiments could indeed be done for QWs as well. 

The ESR spectrum of a LL turns out to be extremely sensitive to the presence and precise 
nature of the spin-orbit coupling mechanism. This is intuitively clear as spin-orbit coupling 
can sometimes destroy spin-charge separation, and nearly always represents the leading term 
breaking SU(2) invariance. Unless this symmetry is broken, the ESR spectrum (including 
electron-electron interactions) is simply a delta peak at the Zeeman energy. Remarkably, 
the spin-orbit coupling in semiconductor QWs will spoil spin-charge separation according to 
refs. [p^pHI, while this does not happen in nanotubes Therefore ESR theory for other 
realizations of the LL (such as nanotubes) is quite different and will be given elsewhere 
Furthermore, the ESR theory of ID spin chains |l5j is also distinct due to the presence of 
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Dzyaloshinskii-Moriya interactions and the absence of charge degrees of freedom. To be spe- 
cific, here we focus on the practically important case of QWs. For the theoretical description 
of the QW including electron-electron interactions, we use the LL model. 

The standard LL Hamiltonian is expressed in terms of charge and spin bosons <p v (y = p, a) 
and their conjugate momenta n„ Q, 

fl=Ey/ dx[K v iil + K-\d,A v f] , (l) 

v=p : a 

with respective effective velocities v v and LL parameters K v . Under Abelian bosonization, 
the right- and left-moving (r = R, L = ±) electron operator for spin a =|, |= ± is 

ip ra (x) = -^=L=exp ^i-y/TT/2[r(j)p + ra<f) a + 9 p + a9 a ]^ , (2) 

where 6 V = — J_ dx'H u and A is a short-distance cutoff- length. To ensure SU(2) invariance 
of eq. (Q), we put K a = 1, while for repulsive interactions, the charge interaction parameter 
K p = g < 1. In addition, to respect Galilean invariance manifest in QWs, the velocities v v 
are taken s&v G — vp and v p — vf/q with the Fermi velocity vf- To simplify notation, below 
we often put vf = Ti = 1. 

The main source of spin-orbit coupling in QWs comes from the Rashba and confinement 
electric fields, and for sufficiently large fields leads to a deformation of each branch of the 
single-particle dispersion relation (l2) . The net result of this deformation is the breaking of 
chiral symmetry, since the Fermi velocities for right- and left-moving electrons with the same 
spin will become different. Therefore, with = vl[ = v± and vri = vl,-\ = «2, we have 
Sv — V2 — Vi > and v\ + v-i = 2vp- A dimensionless measure of the spin-orbit coupling 
strength is then given by A = 5v/2vp <C 1. Concrete estimates for quantum wires yield values 
of the order A ss 0.1 jL3|. In effect, the spin-orbit interaction produces an additional SU(2) 
symmetry-breaking term that can be expressed as 

H' = \Jdx [U p d x (f) a + JIMp] . (3) 

Under an external static magnetic field B along the z axis, we have to add a Zeeman term, 

H z = -B J dxS z (x) = --^LJ dxd x <P„ , (4) 

with QefJ-B = 1 for simplicity. Due to band curvature, the Zeeman term also gives rise to a small 
splitting of the spin-up and spin-down Fermi velocities |jl6| , A<j = \vfi — vfi\/2vf ~ B/2Ef- 
For the most interesting cleaved-edge overgrowth samples, the ID Fermi level is Ep > 20 meV 
Jll| ], and therefore \ a -C A even for strong magnetic fields. Based on this observation, the 
Zeeman-induced splitting of the Fermi velocities is neglected here. 

The ESR intensity at frequency to is then proportional to the Fourier-transformed trans- 
verse spin correlation function, 

I{u)) = j 'dtdxe luJt (S + (t,x)S- (0,0)) , (5) 

where S^(t,x) is the transverse component of the uniform part of the spin operator, = 
Jr + Jl> w ^ tn tne s P m currents Jr^l — V^h /^)i>R,L- The Zeeman term (|]) can now be 



A. De MARTINO and R. EGGER: ESR THEORY FOR INTERACTING ID QUANTUM WIRES 3 



transformed away by shifting the field <j) a and simultaneously II 



/'■ 



Bx 1 B A 

n p^ n p-^7-^2- ( 6 ) 



This is a direct consequence of the fact that the spin-orbit term (||) spoils spin-charge separa- 
tion. Albeit the magnetic field only couples to the spin, it now also affects the charge sector, 
and the combined shift ^ is necessary to transform away the Zeeman term. Of course, the 
fermion operators and the transverse components of the uniform part of spin operator are 
then modified, 

) -> cxp liraBx ^ X ip ra (x) , 



2(1 -A") 



and therefore 



J±(x) -> e* ik ° x J%(x) , J±(x) -> e ±ik ° x j£(x) , u F fc CT = - A 2 ) . (7) 

Since the Hamiltonian remains quadratic in the boson fields, the full problem can be solved 
exactly. We first focus on the T = limit where interaction effects are most pronounced. 
After straightforward but tedious algebra, the fermion propagator and hence the transverse 
spin correlation function are obtained. Using eqs. (g) and (|7|), 



with the same expression but x — > —x for the left-moving part. Here the eigenmode velocities 
iti and U2 are in units of vp: 



Ul - 2 = T2 



g- 2 



1 + 2A 2 T V(5- 2 -l) 2 + 8A 2 ( fl - 2 + l) . (9) 



1/2 



In the absence of spin-orbit coupling, A = 0, these velocities reduce to the spin and charge 
velocities of the LL, u% — * vp and 112 — > Vp/g, while in the non-interacting case, g = 1, we 
recover u\ — > vi and U2 — > t>2- The exponents £^ in eq. (13) are given by 



,±_ t ui ±z 9' 2 + (g' 2 + ^!Zg ± u ^ g' 2 - A2 ) ,■,„■ 

The ESR intensity © follows from eq. (g) and the corresponding left-moving part: 



/(w) = 2 Re / dxdt *, n — TT ( -V' ( ^-tt-V' • (11) 



(27rA) 2 \x — Uit + iAj \x + Uit — ik 

We are unable to compute this double integral in closed analytical form for arbitrary parame- 
ters. However, it has exactly the same structure as the integrals appearing in the computation 
of the spectral function for the spinful LL |l7], |l8). As in that case, one can determine the 
thresholds and the singularities of the ESR spectrum. 

Before that let us briefly analyze some limiting cases. In the non-interacting limit g = 1, 
the exponents ([l0|) simplify to = 1 and = 0. With these values for the exponents the 
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double integral in eq. (^J) can easily be evaluated, and the ESR intensity is uniform over a 
finite range of frequencies with bandwidth Alu = 2BX/(1 — A 2 ) and zero otherwise, 

7H - —\ ■ (12) 

Clearly, in the limit of vanishing spin-orbit coupling (A — > 0), the single delta peak at u> = B 
is recovered. The integrated spectrum is given by 

duj , . IB 1 

Furthermore, in the absence of spin-orbit coupling but with g < 1, it is straightforward to 
recover the delta peak at u — B directly from eq. ( JTT| ) . 

Let us now turn to the general case with g < 1 and A > 0. Following standard argu- 
ments the ESR intensity vanishes identically for to < uik a , with a power-law singularity 
approaching w — u\k a from above, 

I{uj) cx 6(u - u x k a ){u - u 1 k a )^ + ^+^- 1 . (14) 



To see this we change integration variables to s = uit — x and s' = uit + x in eq. (ftl|). In the 
new variables, the integrand exhibits branch cuts in the upper part of the s and s' complex 
plane but is analytic in the lower part. This produces precisely the threshold behavior (14]), 
where the exponent comes from a simple power-counting argument related to a rescaling of s' . 
In order to investigate the behavior near uj = k a u 2 it is more convenient to use the change of 
variables s — u 2 t — x, s 1 = u 2 t + x. In the new variables power counting leads to the following 
power law near oj = k a U2'- 

OC \UJ - U 2 fca|^ +Cl++ ^ _1 • (15) 

However, we cannot prove analytically that the intensity has a second threshold at this fre- 
quency, because the branch cuts of the integrand appear both in the upper and the lower half 
of the s and s' planes. More generally, one can show that because of u 2 > u% it is impossible to 
find a linear transformation, s = axix + a\it and s 1 = a2\x + a 2 2t, such that the position of the 
branch cuts of the integrand produces a Q{k a u 2 — u>) function. However, numerical evaluation 
of e q. (Ill] ) indicates threshold behavior also at the upper edge. Since in general ^ , see 
eq. (10|), the ESR spectrum is markedly asymmetric at both thresholds. The asymmetry is 
quite dramatic, because near to = u\k a , the intensity diverges since the exponent is negative, 
whereas near lo = u 2 k G it vanishes since the exponent is positive. 



The full ESR spectrum can be obtained numerically by direct integration of eq. (11). In 
the numerical integration, a finite cutoff A has to be employed, which leads to several artefacts 
that disappear for A — * 0, namely (i) small oscillations, (ii) finite but very small intensities 
outside the window discussed above, and (iii) negative intensities. Since the time discretization 
should resolve A, numerically one cannot choose arbitrarily small A, and typical results for 
two values of A are shown in fig. [|. Nevertheless, the data in fig. [j] are clearly consistent 
with the general behavior discussed above. This behavior should be compared to both the 
non-interacting case (where the power-law exponents are zero and a constant ESR intensity 
is found in the relevant window) and to the A = case (where one has a (S-peak at the lower 
edge). The interaction-strength-dependence of the exponents close to the lower and upper 
threshold is shown in fig. |2|. The exponents are of different sign, but almost (yet not exactly) 
equal in magnitude. 
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Using ^ ± = £f + £3 and £ = £ + + £ , we find for the integrated ESR spectrum 

7 2^ 7(CJ) ^ 7rr(£+) Ar(£-) + r(i-£+) ( } A 



(16) 



In the non-interacting limit, £~ — ► and £ + — > 2, such that the first term vanishes while the 
second term reproduces eq. (|i~3| ) . Let us finally comment on the effect of thermal fluctuations 
(T > 0). Since the Hamiltonian is quadratic, these can be incorporated by substituting the 
T = propagators in eq. ( |i"l| ) according to 

±iA \ ( ±iA \ 



x^u l t±ikj V sinh — (x =p tiji ± iA) / 

In the non-interacting case, /(w) can then be computed exactly. This solution demonstrates 
that thermal fluctuations wash out the thresholds, broaden the spectrum, and suppress the 
ESR intensity. In the interacting case, these features are also expected. In addition, scaling 
shows that the ESR intensity receives an overall power-law factor T~ 2+ ^. The exponent is 
always positive, vanishes for g = 1, and is generally very small (for A = 0.1, it reaches a 
maximum value of 0.001). 

To conclude, our study of the ESR spectrum resulting from the combined presence of 
interactions and spin-orbit coupling in ID quantum wires reveals dramatic changes from simple 
Fermi liquid modelling. The appearance of asymmetric power-law threshold behavior, with a 
divergence at the lower edge and a vanishing intensity at the upper edge, are manifestations of 
the breakdown of spin-charge separation because of spin-orbit coupling. Such features should 
make an experimental study of our predictions worthwhile and feasible in practice. 
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Fig. 1 - ESR intensity (in arbitrary units) for g = 1/2, A = 0.1, and two values of the cutoff A. The 
threshold at uj = u\k a = 0.988 is indicated by the dashed line, and the numerical result for I[ui) is 
given as the solid curve. The ESR spectrum for A — > terminates at uj = U2k a — 2.011. 
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